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Abstract 

Properties of distributions of the number of trades in different intraday 
time intervals for five stocks traded in MICEX are studied. The dependence 
of the mean number of trades on the capital turnover is analyzed. Correla- 
tion analysis using factorial and H q moments demonstrates the multifractal 
nature of these distributions as well as some peculiar changes in the correla- 
tion pattern. Guided by the analogy with the analysis of particle multiplicity 
distributions in multiparticle production at high energies, an evolution equa- 
tion relating changes in capital turnover and a number of trades is proposed. 
We argue that such equation can describe the observed features of the dis- 
tribution of the number of trades in the stock market. 



1 Introduction 

Statistical analysis of the stock market properties undertaken from the physi- 
cist's point of view has drawn a lot of attention in the last decade. In fact, a 
large fraction of research in the new field of econophysics (UEH^Hl] is devoted 
to questions related to finance. One of the central issues in description of 
the stock price dynamics is understanding the activity pattern characterized 
by the intensity of trading (number of trades in some given time interval). 
In particular, the long-range volatility correlations (one of the most impor- 
tant properties of the financial time series) are directly induced [5] by the 
long-range correlations in trading activity discussed in [0] . Very recently the 
multiscaling of the stock market activity was discussed, in the framework 
of a network approach to complex system dynamics, in [7j. The analysis of 
the properties of the distribution of the number of trades has also allowed to 
prove that the trade-generating process is of distinctly non-markovian nature, 
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so that the time at which the new trade occurs depends on the times of a lot 
of preceding trades [HJ E] . Therefore the study of higher-order correlations is 
of great importance. 

In the present paper we perform a detailed analysis of the distributions of 
the number of trades in a set of chosen intraday time intervals. The technique 
used in our study is borrowed from the analysis of multiparticle production in 
high energy physics, see, e.g., the review papers piU lllj. The analogy drawn 
is between the distribution of the number of trades in some chosen temporal 
interval (finance) and the distribution of the number of produced particles 
in some chosen rapidity (longitudinal momenta) interval, both referred to 
as multiplicity distributions. This technique allows to study correlations in 
the systems of many trades (particles). Both total and genuine (i.e., non- 
reducible to those of smaller groups of trades) correlations are considered for 
sets of any number of trades. The analysis made in the context of multiparti- 
cle production has revealed a number of remarkable properties of multiplicity 
distributions in different rapidity windows, in particular - their multifractal 
nature and a special behavior of genuine correlations governed by a spectac- 
ular dynamics of the particle production process. Therefore it is of obvious 
interest to repeat this analysis in the context of stock price dynamics, which 
is a main issue addressed in the present paper. Its results inspired us to 
formulate a dynamical model of trading activity in financial markets. 

2 Data processing 

The data we used in our analysis is the tick data for five stocks - EESR,MSNG, 
RTKM, SNGS and LKOH, traded at the Moscow International Currency Ex- 
change (MICEX) in the year of 2003. For our study we have chosen five time 
intervals AT: 5 min., 15 min., 45 min. and 495 min. (the trading day at 
MICEX). 

Before turning to the analysis of the multiplicity distributions, let us 
discuss in some more details the above-mentioned analogy between the trade- 
generating process in finance and multiparticle production in high energy 
physics. 

In high energy physics one studies the particle production at some fixed 
collision energy. The higher the energy, the more particles are produced. In 
"financial" terms, one "invests" energy and gets particles from it. In financial 
context we suggest therefore an analogy between the energy and the total 
capital turnover (inflow for buy trades and outflow for the sell ones) in some 
chosen time interval. More precisely, we suggest an analogy between the 
collision energy E co u = ^/s (in the center of mass of colliding particles) and 
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the capital turnover J . The total capital turnover in larger temporal intervals 
is thus disseminated into those in smaller ones corresponding to particular 
distribution of capital in between the trades. 

Let us illustrate the above- described analogy by considering the simplest 
characteristics relating the capital turnover and trade multiplicity distribu- 
tions - a dependence of the mean number of trades in the time interval under 
consideration (ti)at on the capital turnover J~at (in multiparticle dynamics 
the corresponding quantity is thus an energy dependence of mean multiplic- 
ity which has been extensively studied, e.g., in f2lEj)- In Fig- we show 
this dependence in the case of AT = 45 min. 

We see that the mean multiplicity is a concave function of the capital 
turnover J, i.e. (N& T ) oc J a with a < l 1 . For other time intervals the 
dependence is similar. 

Equivalently, for the logarithmic turnover y oc In J~, 



Let us now turn to the analysis of other properties of the trade multiplic- 
ity distributions. The corresponding distributions (normalized at the mean 
number of trades in each interval) for the five stocks considered are shown in 



From Fig. 121 we see, that normalizing the trade multiplicity distributions 
at the corresponding averages brings them to the scale-invariant (i.e. inde- 
pendent of the size of the time window) form. Let us also note that none of 
the standard distributions we have tried (Poisson, negative binomial, gamma, 
inverse gamma) give a satisfactory fit of the data showing that we are dealing 
with some complex trade-generating dynamics. 

The properties of the trade multiplicity distributions are conveniently 
summarized by their suitably defined moments. Among the most important 
characteristics that are widely exploited, in particular, in analyzing the large 
multiplicity events in multiparticle production at high energies, are factorial 
moments (normalized and unnormalized) , cumulants and their ratio, the so- 
called H q moments defined below. 

The unnormalized factorial moments T q with integer ranks q > 1 are 
defined by the following formula: 





Fig. El 



T q = P{n)n(n - l)...(n - q + 1) 



d"G(z) 
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^n multiparticle dynamics this *would correspond to (n) ~ s' 
acteristic to,* e.g., hydrodynamic models of particle production 
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where P{n) is the probability of having n trades (the multiplicity distribu- 
tion) for a given time interval. 

The factorial moments (j2j) are, by defmiton, always positive. It is easily 
seen that the average multiplicity is given by J-\, the dispersion of the distri- 
bution is related to JF 2 etc. The higher is the rank of the moment q, the more 
important is the contribution coming from the tail of the distribution. They 
present therefore an ideal tool for focusing on the properties of the periods 
with intensive trading. 

In equation we have introduced a generating functional for the mul- 
tiplicity distribution 

oo 

G(z) = E^W- (3) 

n=0 

One often considers the normalized factorial moments 

T T 

F = = ^ (A) 

For Poisson distribution (independent trades) all normalized factorial mo- 
ments are equal to 1. Thus a dependence on the overall scale (mean multi- 
plicity) is eliminated. Therefore the moments (HJ) provide a convenient basis 
for extracting information on the structure of correlations. Let us stress, 
that the factorial moments (both normalized and unnormalized) are charac- 
terizing the overall correlation pattern, where the contributions of genuine 
(irreducible) and reducible correlations are mixed. By genuine correlations 
we mean those in which all subsets of the analyzed system are intercon- 
nected. The reducible ones appear in those systems which can be split into 
two or more subsystems not connected to each other. The genuine irreducible 
correlations are described by the cumulants /C: 



d q \nG(z) 
dzi 

For example, if only two-point correlations are present, all cumulants with 
ranks q > 2 are zero, etc. In complete absence of correlations (Poisson) all 
cumulants are zero. 

In what follows we shall focus ourselves on the properties of the normal- 
ized factorial moments F q . In Fig. El we show the dependence of F q on the 
rank q. 

We see, that, apart from the small bending at small q, the logarithms 
of the normalized factorial moments are approximately proportional to the 
rank, In F q oc q, corresponding to an approximate exponential growth of the 
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factorial moments. Let us note in passing, that already this excludes the 
negative binomial distribution, for which the growth of factorial moments 
with q is faster than exponential. 

Of special interest is, of course, a dependence of the moments on the 
temporal resolution AT, revealing the scaling properties of the underlying 
trade multiplicity distributions. In particular, for generic fractal system one 
has \nF q (AT) = a q ln AT. If a q oc q, the dynamics is of monofractal nature. 
More complicated g-dependence of a q indicates that we are dealing with a 
complex multifractal process. The dependence of the normalized factorial 
moments F q on AT is shown, for the moments with ranks in the interval 
q = 1 • • • 10, in Fig. H 

As mentioned above, the scaling pattern is described by the set of the 
slopes in Fig. 0J In Fig we show the ratio of the corresponding average 
slopes to the rank q, i.e., a q /q . 

The simple monofractal scaling would correspond to a q (q)/q = const., 
while from Fig. [3] we clearly see that this is not the case thus pointing out 
to some complex multifractal pattern. In Fig. 0] we observe a characteristic 
strengthening of fluctuations (characterized by the normalized factorial mo- 
ments) with increasing resolution AT. The similar feature is well known in 
quantum chromodynamics [14J 2 . 

As mentioned above, the factorial moments describe the overall (both re- 
ducible and irreducible) correlation pattern. To characterize the irreducible 
content of the multiplicity distribution represented by cumulants, it is con- 
venient to consider the moments H q = K q /T q which represent a share of 
genuine correlations in their overall amount. They are more conveniently 
obtained from the factorial moments (jlj) by the recurrent relation 

9-1 T(a) F F 

The moments H q , plotted in Fig. El show a very interesting behavior: 
they are changing sign at some q ! This phenomenon was first predicted 
theoretically in multiparticle dynamics where it was used as a tool al- 
lowing to uncover the details of the particle generation process. Afterwards 
it was confirmed by experiment [2D E2j- Let us stress, that when working 
with "usual" distributions like Poisson, negative binomial or gamma, the mo- 
ments H q are non-negative. More exactly, they are zero for Poisson, positive 
for negative binomial and non-negative for gamma. To make it possible that 

2 Multifractality of the price-generating process was discussed, in a different context, in 

[Ei ei Ea chi nm 
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H q become negative, one should consider, at least, a mixture of these stan- 
dard distributions, see |2H] ■ Physically, positive cumulants (and thus positive 
H q ) correspond to the effective repulsion between the points, and negative 
ones - to their attraction (clustering of trades). Changing sign means passing 
through the "Poissonian" point, where the corresponding irreducible contri- 
bution is absent. Below we shall consider a phenomenological dynamical 
model of the capital turnover that allows to describe this phenomenon in the 
financial context. 

3 Mathematical model and its predictions 

The above analysis of the stock market data allows to formulate a mathe- 
matical model describing the trade multiplicity distributions inspired by ideas 
originating from high energy physics. As already mentioned, the method of 
correlation analysis discussed below was first suggested in [20] and is widely 
used in the studies of multiparticle dynamics in high energy physics (see, 
e.g., a review article JI]). In multiparticle dynamics one considers multi- 
plicity distributions in some region of phase space at given collision energy. 
Particle creation is described within well-established and well-tested theory 
of strong interactions - quantum chromodynamics. There exist sufficiently 
rigorous equations and well-defined predictions for the characteristics of inter- 
est. Among them are those for the energy dependence of mean multiplicity, 
multifractal properties of the particle production process and specific behav- 
ior of Hg moments. We have observed the similar regularities in the trade 
distributions analyzed above. 

As has been already mentioned, in particular, the change in the sign 
of if^-moments with growing q in multiparticle dynamics was theoretically 
predicted [211] and afterwards confirmed experimentally [2H 122] ■ It has been 
extensively studied theoretically both for the electron-positron collisions [24, 
123 I2E| and for hadron-hadron ones |23j. Let us note, that for standard 
probability distributions such as the negative binomial one this effect does 
not take place. This points out to a quite peculiar underlying dynamics both 
in multiparticle production and trading in the stock market. 

Let us try to write down phenomenological equations for the generating 
functional G(z, y) for the probability distribution of the trade multiplicity at 
given logarithmic turnover y = lnJ' / \ T /J' (j7o= const): 

oo 

G(z,y) = £z"P(n,y) (7) 

n=0 

In Eq. (0) the dependence on y is shown explicitely (while it was not 
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indicated in (j3J)). 

Following the analogy with quantum chromodynamics, we write 



dG 
dy 
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Equation (jHJ) can be interpreted as follows. What we are trying to es- 
tablish is a relation between the distribution of trades in the intervals with 
differing capital turnover. The relation can be described as an evolution 
from smaller to larger capital turnovers. It proceeds through adding "extra" 
turnover, of fractional size x, to the current one with fraction 1 — x forming 
the total turnover given by y - thus the quadratic structure in the RHS of 
(JHJ), describing the inflow to the total turnover. The linear term with neg- 
ative sign corresponds to the outflow from the total turnover and therefore 
it depends on y only. The weight of changing some fraction of the existing 
capital is specified by the kernel K(x) 3 . In equation (jHJ) the constant 70 de- 
scribes the strength of the overall evolution rate in a financial market and 
the integration is over the splitting fraction x. The conservation of turnover 
is due to the arguments In 2 and ln(l — x) in the right hand side. 

Let us assume that the kernel K(x) has the following functional form: 



where c and d are constants. The 1/x term corresponds to the leading role of 
the trades having parametrically small contribution to the capital turnover 4 . 
It is tempting to assume that this effect can be related to the 1/f behavior 
found in the analysis of spectral densities of daily trades [6 . The constant 
and linear terms describe the impact of large trades. In principle, next terms 
of the Laurent expansion of the kernel can contribute as it happens in quan- 
tum chromodynamics. However, in characteristics to be considered, it turns 
out that by keeping the c- and d - terms in the kernel one can conveniently 
parametrize all important effects. The equation (jHJ) with the kernel Q can 
be solved explicitely (see, e.g., [H|). For the mean multiplicity the expo- 
nential dependence on y follows directly if the coupling strength 70 does not 
depend on the capital turnover. Some non-linearity in the exponent can be 
introduced by imposing such a dependence as it happens in quantum chro- 
modynamics (see, e.g., [12] )• The equation (jHJ) also reveals the multifractal 
properties of the distributions [Hj noticed by us in the above analysis. 

3 Note that for technical reasons it is convenient to work with the splitting variable x 
denoting the relative yield of the incoming capital in the newly formed (i.e., total) one. 

4 In particle physics this corresponds to an emission of particles having small energies - 
the so-called "infrared" particles. 




(9) 
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To proceed with iJ g -moments, let us rewrite equation (J7J) as an expansion 



in the normalized factorial moments F q . 



G(z,y) = $»> 9 ^- 

9=0 



(10) 



Assuming that the dependence on the capital turnover y is coming only 
through the average multiplicity (n) 



(n) oc exp(7y), 

we obtain the following equation for the factorial moments 5 : 

i 



dxK(x) 



1 £x" l {l-x)^FiF q . l -F q 
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;i2) 



Equation (|T2*j) is solved iteratively starting from q = 1. At q = 1 one gets 
the following equation for the rate 7 of the growth of mean multiplicity with 
increasing capital turnover (Fq = F\ = 1): 

7 = 7 2 Z" 1 [x 7 + (1 - x) 7 - 1] . (13) 

JO 

From ffTTTf) one can find 7 as a function of 70 and c. With good accuracy 
the corresponding equation can be reduced to an algebraic one: 



7 



7o 2 



C-c 



7 



7T 

IF 



2c- 



7 



where C ~ 0.577 is Euler constant. At 7 1 one has 7 
accurately, 



(14) 
7o- More 



7 ~ 7o 



1 



1 



1 + - (C - c - -d 70 



The magnitude of 7 determines the growth of the average trade multiplicity 
(n) oc exp(7y) with y. In the weak coupling regime 7 is practically equal to 
7o- 

To understand the behavior of the moments H q , let us first consider the 
case of c = d = 0. Then the contribution to the RHS of the evolution 
equation (fT^j) is dominated by small x, and G(y + ln(l — x)) can simply be 
replaced by G(y), so that 



(lnG)' 



7o 2 



1 <ix 



x 



[G(y + lux) - 1] 



7o 2 



dz[G(z) - 1]. (15) 



^Technically the condition of the approximate i^-scaling F' jqF q should hold. 
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Differentiating both sides of (JT5j) over y, we have 



(In G)" = 7o 2 [G-l]. (16) 
Using Eqs (J2J, ©, one gets at 7 ~ 70 

^ = -tl « 4" ( 17 ) 
7 2 g 2 

From Eq. (fTTj) we conclude that the if g -moments are always positive if 
there are small contributions to the capital turnover only. 

Taking into account the remaining terms in the kernel K(x), we get 



H q = ^L_{l + (c-c-\d) q ). (18) 



Now, H q change sign at some value of q. We see, that the zero of the H q 

moments is located at g = (c + \d - Cf \ At larger g the moments H q are 
negative and reach their minimum at some q m in- The larger is the contribu- 
tion of "strong" dealers (large c and (or) d) the lower is the intercept of H q 
with the g-axis. In our analysis of stock data we have just met with such 
a behavior. Let us also note, that from Fig. (jHJ) we see, that the data with 
lower two-point (H 2 ) correlation have a smaller intercept of H q . This im- 
plies that "strong" dealers are more influential in markets (or time intervals) 
characterized by weak interaction with other market participants. Various 
stocks differ in their relative roles - as seen from comparison of various plots 
in Fig. (jnj). It would certainly be of interest to compare these intercepts for 
different markets. This is an important signature of long-range correlations 
involved. 

To make the above analysis more exact, a numerical recursive solution of 
(JI2J) is required. For example, at q = 2 we get the solution in the form 



Jq dxK{x)x 1 {l — x) 1 

2 ~ 2 Jo 1 dxK(x) [x-y + (1 - xfi - 1] - Jo 1 dxK(x) [x^ + (1 - x) 2 ~< - 1] ' 

(19) 

In the above qualitative argumentation we have used, following the anal- 
ogy with QCD, the small - 7 approximation. To get a realistic description of 
financial data in question one has to perform a deeper analysis of the formu- 
lae describing, e.g., the H q moments in the suggested formalism. At present, 
we limit ourselves by these qualitative statements. Quantitative estimate of 
theoretical parameters from the data and its implications will follow. 
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4 Conclusions 



We have analyzed the tick data for five stocks traded at MICEX and found 
some common qualitative features for all of them. The most important ones 
are: 

1. The normalized trade multiplicity distributions are, to a good accuracy, 
scale-invariant with respect to the size of the time window considered. 

2. The average multiplicity of trades is a concave function of the total 
capital turnover. 

3. The multiplicity distributions demonstrate a multifractal behavior. 

4. The cumulant (and therefore H q ) moments change their sign at some 
rank. 

These features are strongly reminiscent of the properties of particle multi- 
plicity distributions studied in high energy physics where the well-developed 
theory of quantum chromodynamics describes (and moreover predicts) them 
quite successfully. Using this analogy, we have proposed the phenomenologi- 
cal model of trades multiplicity distribution with the equation similar to that 
arising in quantum chromodynamics. This equation provides a qualitative 
description of the above-mentioned properties. Its interpretation in terms of 
inflow and outflow of the capital turnover, providing some insight into the 
dynamics of financial markets, is considered. 

The work was supported by RFBR grants 04-02-16880, 04-02-16445-a, 
02-02-16779 and the Scientific Schools Support Grant 1976.2003.02. 
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Figure 3: Normalized factorial moments F q ; AT = 5 min (red), 15 min 
(blue); 45 min (brown); 495 min (black). 
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